All-sky search for gravitational-wave bursts in the second joint LIGO-Virgo run 
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We present results from a search for gravitational-wave bursts in the data collected by the LIGO 
and Virgo detectors between July 7, 2009 and October 20, 2010: data are analyzed when at least 
two of the three LIGO- Virgo detectors are in coincident operation, with a total observation time of 
207 days- The analysis searches for transients of duration < 1 s over the frequency band 64-5000 Hz, 
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without other assumptions on the signal waveform, polarization, direction or occurrence time. All 
identified events are consistent with the expected accidental backgrotmd. We set frequentist upper 
limits on the rate of gravitational-wave bursts by combining this search with the previous LIGO- 
Virgo search on the data collected between November 2005 and October 2007. The upper limit on 
the rate of strong gravitational- wave bursts at the Earth is 1.3 events per year at 90% confidence. 
We also present upper limits on source rate density per year and Mpc 3 for sample populations of 
standard-candle sources. As in the previous joint run, typical sensitivities of the search in terms of 
the root-sum-squared strain amplitude for these waveforms lie in the range ^5x 10“ 22 Hz” 1 / 2 to 
~ 1 x 10“ 20 Hz -1 / 2 . The combination of the two joint nans entails the most sensitive all-sky search 
for generic gravitational-wave bursts and synthesizes the results achieved by the initial generation 
of interferometric detectors. 

PACS numbers: 04.80.Nn, 07.05.Kf, 95.30.Sf, 95.85.Sz 


I. INTRODUCTION 

Astrophysical sources of transient gravitational waves 
(duration of < 1 s) [1] include merging compact bi- 
nary systems consisting of black holes and/or neutron 
stars [2, 3], core-collapse supernovae [4], neutron star 
collapse to black holes [5], star-quakes associated with 
magnetar flares [6] or pulsar glitches [7], cosmic string 
cusps [8], and other violent events in the Universe. Since 
many classes of gravitational- wave (GW) bursts cannot 
be modeled well - if at all - a search for those sources 
must be sensitive to the widest possible variety of wave- 
forms, 

This paper reports on a search for GW bursts occurring 
during the second joint observation run of the LIGO [9] 
and Virgo [10] detectors, which took place in 2009-2010. 
This search makes no prior assumptions on source sky lo- 
cation, signal arrival time, or the waveform itself. Event 
rate upper limits from long-term searches of this category 
have been derived with networks of resonant bar detec- 
tors with spectral sensitivity limited to around 900 Hz in 
1997-2000 [11, 12] and in 2005-2007 [13, 14]. Networks 
of interferometric detectors set more stringent upper lim- 
its for GW bursts on a wider bandwidth using the LIGO 
detectors in 2005-2006 [15-17] and during the first joint 
observation of LIGO and Virgo detectors in 2007 [18]. 

This second joint LIGO-Virgo search for GW bursts 
analyzed the frequency band spanning 64-5000 Hz. We 
achieved a frequency-dependent sensitivity comparable 
to or better than that of the first joint run, and accu- 
mulated ~ 207 days of observation time interlaced with 
periods of installing or commissioning major hardware 
upgrades. Moreover, for the first time a low-latencv 
analysis was run with the goal of providing triggers for 
electromagnetic follow-ups of candidates by robotic opti- 
cal telescopes [19], radio telescopes, and the Swift satel- 
lite [20, 21]. In this paper we focus on the final results of 
the GW stand-alone search, which found no evidence for 
GW bursts. 

This paper is organized as follows. In Section II we 
describe the second joint scientific run: we report on the 
LIGO and Virgo instrumental upgrades with respect to 
the first run and on data quality studies. In Section III 
we give a brief overview of the search: the search algo- 


rithm, background estimation, the simulations and the 
calibration uncertainties. The signal models (GW wave- 
forms and source populations) we tested are described in 
Section IIIC. The results of the search are presented in 
Section IV, and astrophysical implications are discussed 
in Section V. The Appendices provide additional details 
on data characterization and analysis methods. 


II. SECOND LIGO-VIRGO SCIENCE RUN 

The network of detectors used in this search comprises 
the two LIGO 4 km interferometers, denoted ‘'HI” (lo- 
cated in Hanford, WA) and ‘LI” (Livingston, LA), as 
well as the Virgo 3 km interferometer, denoted ‘VI” 
(close to Pisa, Italy) l . 

The LIGO detectors operated from July 7, 2009 to 
October 20, 2010 in their sixth science run (S6). The 
Virgo detector operated from July 7, 2009 to January 8, 
2010 in its second science run (VSR2) and again from 
August 11 to October 20, 2010 in its third science run 
(VSR3). 

As in the first joint LIGO-Virgo run [18, 22], the oper- 
ation of three differently oriented and widely separated 
detectors allows for reasonably complete coverage of the 
sky for at least one gravitational-wave polarization com- 
ponent as well as the recovery of some source character- 
istics such as sky location [19, 23, 24]. 


A. Detector Upgrades 

Before the beginning of the runs, several detector hard- 
ware upgrades were implemented in order to prototype 
new subsystems planned for the next generation of de- 
tectors, referred to as “advanced detectors” [25, 26], ex- 
pected to start observations in 2015. The upgrades of 


1 The 2 km detector at the Hanford site (H2) was decommissioned 
before the second joint LIGO-Virgo run. During previous runs, 
the latter detector was mainly used to enforce additional event 
selection criteria by taking advantage of the special relationship 
for GW signals from the co- located interferometers HI and H2. 
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FIG. 1: Noise spectra for the three LSC- Virgo detectors 
achieved during S6-VSR2/3. 


the LIGO detectors for S6 include a higher power 35 W 
laser, the implementation of a DC readout system, a new 
output mode cleaner, and an advanced LIGO seismic iso- 
lation table [27]. The upgrades of the Virgo detectors 
were achieved in two steps. For VSR2, Virgo operated 
with a more powerful laser and a thermal compensation 
system. Virgo then went offline to undertake a major up- 
grade of its suspensions, mainly in the installation of new 
test masses consisting of mirrors hung from fused silica 
fibers [28]. Virgo resumed observations in August 2010 
with VSR3. Best sensitivities, in terms of noise spec- 
tral densities, of the LIGO and Virgo detectors achieved 
during their second joint run (henceforth defined as S6- 
VSR2/3), as a function of signal frequency, are shown in 
Figure 1. 


B. Data Quality 

To mitigate the consequences of new hardware instal- 
lations and detector commissioning during this run, sig- 
nificant effort has been made to identify and charac- 
terize instrumental or data acquisition artifacts, periods 
of degraded sensitivity, or an excessive rate of transient 
noise due to environmental conditions [29]. During such 
times, the data were tagged with Data Quality Flags 
(DQFs). Following the same approach used in previous 
searches [16-18], these DQFs are divided into three cate- 
gories depending on their impact on the search and on the 
understanding of the behavior of the detector. A further 
description of DQF categories is presented in Appendix 
A. 

After DQFs have been applied, the total analyzable 
time for the S6-VSR2/3 run is 242.8 days for HI, 220.2 
days for LI, and 187.8 days for VI. 


III. SEARCH OVERVIEW 

In this analysis, we considered all four available detec- 
tor network configurations: the three detector network, 
HI LI VI, and the three combinations of detector pairs, 
H1L1, HI VI and LI VI. We decided a priori to search for 
GW bursts in the entire available time of three-fold ob- 
servation and in the remaining exclusive times of the two- 
fold networks. Table I reports the total (non-overlapping) 
coincident observation time for each configuration of de- 
tectors searched for GW signals. Information about dis- 
tinct sub-periods of the run may be found in Appendix C. 

Due to the commissioning breaks and installation ac- 
tivities described in Section II, the total observation time 
is dominated by 2- fold configurations. 


network 

H1L1V1 

H1L1 

L1V1 

H1V1 

total 

observation time [days] 

52.2 

84.5 

28.9 

41.0 

206.6 


TABLE I: Mutually exclusive observation time for each de- 
tector configuration after the application of category 2 DQFs 
(see Appendix A for the definition of data quality flags and 
their categories). 

The useful frequency band is limited to 64-5000 Hz by 
the sensitivity of the detectors and by the valid range of 
data calibration. For computational reasons, the event 
search was performed separately in two suitable bands, 
64-2048 Hz and 1600-5000 Hz, overlapping to preserve 
sensitivity to events with spectral power at intermediate 
frequencies. The analysis of the events (including the 
tuning of the search) was performed independently on 
each configuration of detectors and on three sub-bands, 
namely 64-200 Hz, 200-1600 Hz and 1600-5000 Hz, by 
classifying the found events according to their recon- 
structed central frequency. The motivation for this band 
splitting is to tune the search within event sets of homo- 
geneous glitch behaviour. 

A. Search Algorithm 

This search is based on the coherent WaveBurst (cWB) 
algorithm [30], which has been used since LIGO’s fourth 
science run in various searches for transient GWs [16- 
18, 31]. 

The cWB analysis is performed in several steps. First, 
detector data is decomposed into a time-frequency repre- 
sentation and then whitened and conditioned to remove 
narrow-band noise features. Events are identified by 
clustering time-frequency pixels with significant energy 
which is coherent among detectors and characterized us- 
ing test statistics derived from the likelihood (which is 
also a measure of the signal energy detected in the net- 
work and is calculated as described in [16]). The primary 
statistics are the network correlation coefficient ce, which 
is a measure of the degree of correlation between the de- 
tectors, and the coherent network amplitude r/, which is 
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proportional to the signal-to-noise ratio and is used to 
rank events within a homogeneous sub-period. 

Both of these statistics are described in detail in [16]. 
The application of the event selection criteria is thor- 
oughly described in [18, 32]. 

Any gravitational-wave candidate event detected by 
cWB is subject to additional data-quality vetoes based on 
statistical correlations between the GW data channel and 
environmental and instrumental auxiliary channels; a sig- 
nificant correlation indicates the event may have been 
produced by local noise. Further details can be found in 
Appendix B . 


B. Background Estimation and Search Tuning 

A sample of “off-source” (background) events is re- 
quired to determine the selection thresholds to re- 
ject noise-induced events contaminating the “on-source” 
(foreground) measurement. We estimate the distribution 
of background events by performing the analysis on time- 
shifted data, typically in ~ 1 s steps. The shifts minimize 
the chance of drawing an actual GW into the background 
sample. To accumulate a sufficient sampling, this shifting 
procedure is performed hundreds or thousands of times 
without repeating the same relative time shifts among de- 
tectors. Background events corresponding to times which 
are flagged by data quality studies are discarded, just as 
an event candidate from the foreground would be. 

Due to the different characteristics of the background 
noise for the various sub-periods between commissioning 
breaks and for the different frequency bands and net- 
works, the thresholds on cc and 77 are tuned separately 
for each homogeneous sub-period. Moreover, we consider 
the action of conditional DQFs (Category 3 DQFs, see 
Appendix A) on the event significance, by introducing 
a new ranking scheme which assigns lower significance 
to events flagged by such DQFs. More details on this 
procedure are reported in Appendix D. 

The thresholds reported in Table VII are selected to 
require a false alarm rate (FAR) ~ 1/(8 yr) 0 r lower per 
frequency band (to be considered as a warning thresh- 
old), as described in Appendix D. By considering the 
union of all searches performed (network configurations, 
sub-periods, and frequency bands) we estimate an overall 
false alarm probability (FAP) of ~ 15%. 

The method to measure the statistical significance 
of any “on-source” candidate GW event which passes 
the aforementioned false alarm thresholds was decided 
a priori , namely performing additional independent time 
slides to increase the statistics of the background esti- 
mates. 


C. Simulated Signals and Detection Efficiencies 

In order to test the sensitivity of our search to 
gravitational-wave bursts, we add (“inject”) various ad- 


hoc software signals, both polarized and un-polarized, 
to the detector data and measure the detection effi- 
ciencies of the search. The injected waveforms can be 
parametrized as: 


M*)" 

= Ax 

retail 

2 

X 

H + (tj 

h x (t) 


Ol 


H x (t) 


where A is the amplitude, a the ellipticity 2 and H + / x 
are the waveforms for the two independent polarizations. 
In this search, we investigated elliptically polarized sig- 
nals (i.e. ol uniformly chosen in [0,1]), as well as sets of 
only linearly or circularly polarized waves (or fixed to 0 
or 1, respectively). A variety of GW signal morphologies 
spanning a wide range of signal durations, frequencies 
and amplitudes were tested. See Figure 2 for a sample of 
representative waveforms from various families and Ta- 
bles II, III, IV for the chosen waveform parameters. 

The injected waveform families include: 

• Sine- Gaussian: 

H+(t) = exp {~t 2 /r 2 ) sin(27r f Q t) (3.2) 

H x {t) = exp (— £ 2 /t 2 ) cos(27t/o£) (3.3) 

where r = Q/(\/27r/ 0 ). We consider waveforms 
of this type with central frequencies /o chosen be- 
tween 70 to 5000 Hz and quality factors Q = 
3,9, 100. Sine-Gaussian waveforms with a few cy- 
cles are qualitatively similar to signals produced by 
the mergers of two black holes [2]. 

• Gaussian: 

H + (t)=exp{-t 2 /r 2 ) (3.4) 

H x (t)= 0 (3.5) 

where the duration parameter r is chosen to be one 
of 0.1, 1.0, 2.5, or 4.0 ms. 

• Ring- down wavefor,ns: 

H+(t) = exp (-t/r) sin(27 rf 0 t) (3.6) 

H x ( t ) = exp (-t/r) cos(27r f 0 t) (3.7) 

We use several central frequencies from 1590 Hz to 
3067 Hz, and decay times r = 0.2 s or Q = 9. Ring- 
downs can occur in the end stages of black hole 
binary mergers. Longer duration ring-downs are 
also similar to signals predicted from the excitation 
of fundamental modes in neutron stars [33], 


2 For binary sources, the ellipticity i« the cosine of the source in- 
clination angle, i.e. the angle between the source rotational axis 
and the line of sight to Earth. 
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FIG. 2; Representative waveforms injected into data for simulation studies. The top row is the time domain and the bottom 
tow is a time-frequency domain representation of the waveform. From left to right: a 361 IIz Q = 9 sine-Gaussian, a r = 4.0 ms 
Gaussian waveform, a white noise burst with a bandwidth of 1000-2000 Hz and characteristic duration of r = 20 ms and, 
finally, a ringdown waveform with a frequency of 2000 Hz and r = 1 ms. 


• Band-limited ■ white noise signals : 

The polarization components are bursts of uncorre- 
lated band-limited white noise, time shaped with a 
Gaussian profile; H+ and H x have — on average — 
equal RMS amplitudes and symmetric shape about 
the central frequency (see Figure 2). 

• Neutron star collapse waveforms : 

For a comparison with previous searches [17, 18], 
we considered numerical simulations by Baiotti et 
al. [5], who modeled neutron star gravitational col- 
lapse to a black hole and the subsequent ring-down. 
As in previous searches, we chose the models D1 (a 
nearly spherical 1.67 M 0 neutron star) and D4 (a 
1.86 M 0 neutron star that is maximally deformed 
at the time of its collapse into a black hole) to rep- 
resent the extremes of the parameter space in mass 
and spin considered in the aforementioned work. 
Both waveforms are linearly polarized (H x = 0) 
and their emission is peaked at a few kHz. 

The simulated signals were injected with many ampli- 
tude scale factors to trace out the detection efficiency as 
a function of signal strength. The amplitude of the sig- 
nal is expressed in terms of the root-sum-square strain 
amplitude ( h rss ) arriving at the Earth, defined as: 


1 J IM0l 2 + IM0l 2 <ft (3.8) 

The signal amplitude at a detector is modulated by the 
detector antenna pattern functions, expressed as follows: 

Adet(t) = F + (e,4>,^)/i+(t) + F*(0, (3.9) 

where F+ and F x are the antenna pattern functions, 
which depend on the orientation of the wavefront relative 
to the detector, denoted here in terms of the sky position 


(0,$), and on the polarization angle %/k The sky po- 
sitions of simulated signals are distributed isotropically 
and polarization angles are chosen to be uniformly dis- 
tributed. The detection efficiency is defined as the frac- 
tion of signals successfully recovered using the same se- 
lection thresholds and DQFs as in the actual search. The 
detection efficiency of the search depends on the network 
configuration and the selection cuts used in the analysis. 

Detection efficiencies for the H1L1V1 network for se- 
lected waveforms as a function of signal amplitude h TSS 
and as a function of distance (for the D1 and D4 wave- 
forms from Baiotti et al. [5]) are reported in Figures 3 
and 4 , respectively. As in the previous joint run, typical 
sensitivities for this network in terms of h Tsa for the se- 
lected waveforms lie in the range ^5x 10 -22 Hz” 1 / 2 to 
~lx 1(T 20 Hz -1 / 2 ; typical distances at 50% detection 
efficiency for neutron star collapse waveforms lie in the 
range ~ 50 pc to ~ 200 pc. 

Two convenient characterizations of the sensitivity, 
the /i rea at 50% and 90% detection efficiency ( h jjj° % and 
h^° respectively) are obtained from fitting the efficiency 
curves and are reported in Tables II, III, and IV for the 
various families. Notice that the 3-fold network, H1L1V1, 
has a better sensitivity than the weighted average over 
all networks: 2-fold networks have ~ 3/4 of the analyzed 
live time, but feature a lower sensitivity. 


D. Systematic Uncertainties 

The most relevant systematic uncertainty in the astro- 
physical interpretation of our results is due to the cal- 
ibration error on the strain data produced by each de- 
tector [34, 35]. The effect of calibration systematic^ on 
network detection efficiency has been estimated by dedi- 
cated simulations of GW signals in which the signal am- 
plitude and phase at each detector is randomly jittered 


9 


f 0 [Hz] 

— 235 L 
■ •«> 235 E 

— 1304 L 

— 1 304 E 

10 22 10' 21 10 20 

hrss, [strain/^ Hz] 




fo [Hz] 

— 2090 L 
* 2090 C 

2590 L 

— 2590 C 

10“ 10‘ 21 lO 20 

hrss, [strairA Hz] 




c 

0 

1 


5 

Q 



Af t 

Hz H2 ms 
100 ICO 100 

a 1000 10 100 
1000 1000 10 
- * 200C 1000 10 
3500 100 100 
3500 1000 10 


10‘ 21 10' 20 
hrss, [si rain/s Hz] 


FIG. 3: Detection efficiency for selected waveforms as a 

function of signal amplitude /i fS s for the HI LI VI network. 
Top: Comparison of detection efficiency for linear (L) and 
elliptical (E) sine-Gaussians with central frequencies of 235 
and 1304 Hz. Middle: Comparison of detection efficiency for 
linear (L) and circular (C) ring-down signals with frequencies 
of 2090 and 2590 Hz. Bottom: Detection efficiency for white 
noise bursts with frequency spanning between 100 and 4500 
Hz. 


FIG. 4: Efficiency for the H1L1V1 network as a function of 
distance for the D1 and D4 waveforms predicted by polytropic 
general-relativistic models of neutron star collapse. 


be parametrized as an overall shift of the detection effi- 
ciency curves along the signal strength axis. The largest 
effect over the injected signal waveforms was a 8% in- 
crease of the fr ra8 amplitude at fixed detection efficiency 
3 . To produce the astrophysical limits shown in Section 
IV, we use the reduced detection efficiency curves ob- 
tained by shifting the original fits from Subsection III C 
and the results in Tables II, III, and IV to 8% larger h rsf3 
values. 


IV. SEARCH RESULTS 

The on-source data have been analyzed following the 
procedures tuned through the investigation of the off- 
source background sample, as described in Appendix D. 
No on-source event has been found above the threshold 
false alarm rate of once in 8 years per frequency band, 
and the distribution of on-source events is in agreement 
with the measured background. Table V lists the five 
most significant on-source events, as ranked by their In- 
verse False Alarm Rate (IFAR = 1/FAR), and taking 
into account the trial factor due to the three indepen- 
dent searches performed on the disjoint frequency bands. 

In addition to the events reported in Table V, this 
search also detected an on-source event showing a chirp- 
ing waveform compatible with a compact binary coales- 
cence at a signal-to- noise ratio ~ 17 in the H1L1V1 net- 
work. This event was first identified by a low-latency 
burst search within minutes of its occurrence on Septem- 


according to the modeled distribution of calibration er- 
rors for that detector. The resulting network detection 
efficiency marginalizes the effect of the systematic uncer- 
tainties over the observation time. The main effect can 


3 Note that, due to an incomplete knowledge of the actuation res- 
onances in [3000,4000] Hz band of the Hanford detector, very 
conservative assumptions have been made on calibration uncer- 
tainties; the networks including HI in that frequency band fea- 
ture a large efficiency loss due to calibration systematics of 24%. 
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2477 

100 

34.6 


46.3 

233.5 



2756 

9 

44.2 


56.8 

389.8 



3067 

3 

74.1 


81.7 

600.0 



3067 

9 

64.6 


78.0 

499.6 



3067 

100 

41.1 


53.8 

278.2 



3413 

9 

65.7 


80.0 

510.4 



3799 

9 

81.7 


99.3 

719.9 




TABLE II: Values of and h™, % (for 50% and 90% de- 
tection efficiency at the chosen thresholds of 1/(8 yr) per fre- 
quency band), in units of 10 _22 Hz“ 1,2 , for linear and ellip- 
tical sine-Gaussian waveforms with the central frequency f Q 
and quality factor Q. The center two columns are the 
for linear and elliptical waveforms during the total S6 period 
measured for the HI LI VI network. The rightmost columns 
report the values of and /i? s ° 3 % over the whole S6-VSR2/3 
for the combined results (i.e. averaged over time) from all the 
networks. 


ber 16, 2010 and was thoroughly investigated in follow-up 
studies* Its Inverse False Alarm Rate was estimated at 
1.1 yr from comparison with the burst reference back- 
ground over all trials. After the completion of the anal- 
ysis, this event was revealed to be a blind hardware in- 
jection [36] intended as an end-to-end test of the search 





all networks 


/ 

r 

Linear 

Circular 

[Hz] 

M 

as* 

AS? 

l50% 
fL rs a 

i,90% 

N'rss 

2000 

1.0 

47.3 

288 

34.8 

78.9 

2090 

200 

42.9 

218 

31.7 

66.0 

2590 

200 

52.2 

255 

39.1 

79.5 

3067 

0.65 

91.9 

546 

72.9 

569 


TABLE III: Values of tig* and h 9 ™ (for 50% and 90% de. 
tection efficiency at the chosen thresholds of 1/(8 yr) per fre- 
quency band), in units of 10“ 22 Hz~ 1//2 , for linearly and cir- 
cularly polarized ring-downs characterized by parameters / 
and r. 


/low A / r 

[Hz] [Hz] [ms] 

H1L1V1 

1,50% 

'frrss 

all networks 

l50% l90% 

'^rss '*tss 

100 100 100 

8.1 

11.5 91.2 

250 100 100 

7.5 

10.5 43.1 

1000 10 100 

15.5 

22.5 93.6 

1000 1000 10 : 

30.5 

39.7 130 

1000 1000 100 

76.8 

76.7 492 

2000 100 100 

35.7 

40.3 193 

2000 1000 10 

55.6 

63.1 211 

3500 100 100 

71.8 

90.3 332 

3500 1000 10 

114 

125 371 


TABLE IV: Values of and h? s ° s % (for 50% and 90% de- 
tection efficiency at the chosen thresholds of 1/(8 yr) per 
frequency band), in units of 10“ 22 Hz -1/2 , for band-limited 
white noise waveforms characterized by parameters /i ow , A/, 
and r. 


for transient signals 4 . As such, the event was removed 
from the final results. 


4 Signal injections were performed via direct excitation of the 
interferometer mirror test masses. Some of these hardware in- 
jections were intended to mimic a coherent GW excitation across 
the network and to provide an end-to-end verification of the de- 
tector instrumentation, the data acquisition system and the data 
analysis software. In addition to those, a blind injection chal- 
lenge was realized consisting of injecting a few simulated signals 
at times not announced to the collaborations. This was done 
for the purpose of testing the data analysis pipelines and event 
validation protocols. 



11 


IFAR [yr] 

freq. band 

network 

SNR 

FAP 

0.64 

0.2-1.6 kHz 

H1L1 

11 

0.59 

0.36 

64-200 Hz 

H1L1V1 

19 

0.47 

0.28 

0.2-1. 6 kHz 

H1L1 

12 

0.33 

0.19 

0.2-1. 6 kHz 

H1L1 

10 

0.35 

0.17 

1.6-5 kHz 

HI VI 

9 

0.24 


TABLE V: The five most significant events present in the on- 
source data. IFAR is the Inverse False Alarm Rate [yrj of the 
event in the entire search, SNR is the signal-to-noise ratio in 
the whole network, and FAP is the false alarm probability 
(probability of getting at least as many accidental events as 
those observed with IFAR > the value reported in the first 
column). 


A. Upper Limits 

The new null result can be combined with the previous 
ones from the latest scientific runs by LIGO and Virgo 
[16-18] to complete the results achieved by initial gener- 
ation interferometric detectors. 

Assuming a Poisson distribution of astrophysical 
sources and in the special case of no sur\ iving candidate 
events, the 90% confidence upper limit is computed as 

iWo = Y" 3 t > (4.1) 

Z-rfc 

where 2.3 — - log(l - 0.9), th and T* are respectively the 
detection efficiency and the observation time of the differ- 
ent network configurations in homogeneous sub-periods 
of observation k [37], including all available LIGO and 
LIGO-Virgo observations since November 2005 [16-18]. 

Figure 5 shows the upper limits on the rate of 
gravitational-wave bursts at the Earth as a function of 
signal strength (/i rss ) for selected sine-Gaussian wave- 
forms. The second joint LIGO-Virgo run increases the 
previous total observation time by roughly 50%, total- 
ing 1.74 yr. Therefore, the resulting 90% upper limit 
on the rate for strong signals (asymptotic behaviour for 
6^ -» 100%) decreases from 2.0 to 1.3 yr -1 for the 64 - 
1600 Hz band (from 2.2 to 1.4 yr -1 for the band above 
1.6 kHz). 

The results can also be interpreted as limits on the 
rate density of GW bursts (number per year and per 
Mpc 3 ) assuming a standard-candle source, isotropically 
distributed, as previously reported in [18]. Denoting by 
/iq the average value of the GW squared amplitude h* SH 
at a fiducial distance ro from the source, the energy con- 
verted to GWs is 

#gw = -Q- r o fo hi . (4.2) 

where fo is the central frequency of GW emission. 

Considering a population of standard- candle sources 
randomly oriented with respect to the Earth and at a 



70 Hz 
100 Hz 
235 Hz 
361 Hz 
1053 Hz 
2000 Hz 
3067 Hz 


FIG. 5: Upper limits at 90% confidence on the rate of 
gravitational-wave bursts at Earth as a function of h TaB signal 
amplitude for selected sine-Gaussian waveforms with Q = 9. 
The results include all the LIGO and LIGO-Virgo observa- 
tions since November 2005. 


distance r 0 , we can interpret the h q as the average GW 
squared amplitude impinging on the Earth (e.g. averaged 
over source parameters such as inclination angle) . Equar 
tion 4.2 can then be used to estimate /iq(^gw* /o> r o) and 
in particular sets the inverse proportionality between the 
average h TSS at Earth and source distance r: hr = horo. 
Assuming a uniform distribution in the sky and in time 
of these standard-candle sources, the expected rate of 
detections is 

Ndet - 47t7^T [ drr 2 e(r) 

Jo 

= An7ZT(h 0 r 0 ) 3 [ dh fc“ 4 c (fc). (4.3) 
Jo 

where K is the rate density of the standard-candle 
sources, T the overall observation time, and e(h) the de- 
tection efficiency as measured by our simulations. 

Hence, the 90% confidence upper limit on rate density 
71 of such standard-candle sources is 

2.3 

90% 47rT’(/io r o) 3 / 0 °°d/i/i^ 4 e(/i) 

The resulting 71qq% is dominated by the part of the de- 
tection efficiency curve at small GW amplitude h. Due 
to the relative orientations of the LIGO-Virgo detectors, 
detection efficiency curves for linearly polarized sine- 
Gaussian waveforms are approximately the same of those 
for elliptically polarized ones; the numerical values of 
7lg 0 % are close within a few percent for both source mod- 
els. 

Figure 6 shows the rate density upper limits of 
sources as a function of frequency. This result 
can be interpreted in the following way: given a 

standard-candle source with a characteristic frequency 
/ and energy Eq w, the corresponding rate limit is 

7W/)(M e c 2 /£ GW ) 3/2 yr-'Mpc" 3 
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Frequency [Hz] 


FIG. 6: Rate limit per unit volume for standard-candle 
sources at the 90% confidence level for a linearly polarized 
sine-Ga,ussian standard-candle with Equ = M®c 2 . Within 
an accuracy of a few percent, the same numerical results 
hold also for sources emitting circularly polarized GWs, which 
would subsequently appear elliptically polarized at the Earth. 
In this Figure, all LIGO and LIGO- Virgo observations since 
November 2005 have been combined together. 



Linear Q=9 
sine-Gaussian 

Elliptical 0=9 
sine-Gaussian 

Linear 

Ring-down 

Circular 

Ring-down 

Band-Jimited 
white noise 


V. SUMMARY AND DISCUSSION 


This paper reports the results achieved by the LIGO 
and Virgo detectors in the search for GW transients of 
duration < 1 s, without assumptions on the signal wave- 
form, polarization, direction or arrival time. 

Three detectors were operating at the Hanford, Liv- 
ingston and Pisa sites during the second joint observa- 
tion of LIGO and Virgo in 2009-2010. The detectors im- 
plemented hardware upgrades in order to prototype new 
subsystems planned for the upcoming advanced detec- 
tors. The resulting sensitivities to GWs were comparable 
to those achieved during the first LIGO-Virgo rim. The 
main contribution of the second run is a 50% increase in 
accumulated observation time. 

No event candidates were found in this search. We 
set better upper limits on the rate of gravitational-wave 
bursts at Earth and on the rate density of burst sources 
per unit time and volume. These limits combine all avail- 
able information from the LIGO-Virgo joint runs and set 
the state-of-the-art on all-sky searches for transient grav- 
itational waves of short duration. 

The reported h TSS amplitude of the GW at Earth can 
be converted into the energy' emitted by a source at some 
fiducial distance assuming a simple model as in Equa- 
tion 4.2. For example, the energy emitted in gravita- 
tional waves in units of solar masses at a distance of 10 
kpc and considering measured h raa at 50% detection ef- 
ficiency (Table II) is ~ 2.2 * 10“ 8 M G for signal frequen- 
cies near 150 Hz (5.6 • 10~ 2 M o at 16 Mpc). These GW 
energies, though obviously depending on the signal fre- 
quency, are approximately constant over different polar- 
ization models of the GW emission, including linearly 
polarized sources, circularly polarized sources and un- 
polarized emission with random polarization amplitudes 
(see Tables II, III, and IV). 


FIG. 7: Typical GW energy in solar masses at 50% detection 
efficiency for standard-candle sources emitting at 10 kpc for 
the waveforms listed in Tables II, HI, and IV considering the 
H1L1V1 network and the LIGO-Virgo observations since July 
2009. 


The typical GW energy in units of solar masses for 
LIGO-Virgo observation is shown in Figure 7 computed 
with Equation 4.2 using the measured h r8 s at 50% detec- 
tion efficiency for the tested waveforms assuming a stan- 
dard candle source emitting at a distance of 10 kpc. The 
mass scales with the square of the fiducial distance and 
the results are robust over the very wide class of wave- 
forms tested. As expected, the GW energy is strongly de- 
pendent on the spectral sensitivity of the network, with 
a negligible dependence on the specific waveform charac- 
teristics. 


The long baseline interferometric detectors LIGO and 
Virgo are currently being upgraded to their advanced 
configurations, and the next joint observation is planned 
for 2015. Another advanced detector, LCGT [38, 39], is 
being built in Japan, and there are proposals to realize 
an additional advanced LIGO detector outside the USA. 
These advanced detectors should achieve strain sensi- 
tivities a factor of ten better than the first-generation 
detectors. For example, at design sensitivity these de- 
tectors should detect a typical core-collapse supernova 
anywhere in the galaxy [40] and will be able to put con- 
straints on extreme scenarios for core collapse supernovae 
within the Local Group [4, 41]. Other possible short du- 
ration sources, such as the merger of very high mass stel- 
lar black hole binaries, could be visible at distances ex- 
ceeding 1 Gpc, During advanced detector observations, 
gravitational-wave detections are predicted to occur on 
a regular basis [42], thus greatly expanding the field of 
gravitational-wave astrophysics. 
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Appendix A: Data Quality Flags 

Data Quality Flags (DQFs) are intended to indicate 
periods of data taking which suffer from environmental 
and instrumental effects inducing noise into the data [29]. 
We followed the DQF strategy used in previous searches 
[16-18], organizing DQFs into 3 categories. The differ- 
ent categories reflect the level of understanding of the de- 
tectors’ performances as well as of the relation between 
disturbances in the data set and environmental or instru- 
mental causes. 

Category 1 DQFs mark segments of time (typically 
more than tens of seconds) when disturbances make anal- 
ysis unfeasible. Data segments remaining after their ap- 
plication are used in the analysis. 

Category 2 DQFs axe connected to well-understood 
short duration (typically a few seconds) periods of noise 
transients. Data segments flagged by this category can be 
used for data conditioning and noise property estimation, 
but events emerging from these periods are discarded as 
very likely originating from the detector environment. 

Finally, Category 3 DQFs denote periods that are 
only v’eakly correlated to environmental and instrumen- 
tal monitors. Such cuts are not reliable enough to be used 
as unconditional cuts. When applied to events generated 
by the search algorithm, they would reject a significant 


fraction (in extreme cases up to 15-20%) of data. Their 
use is limited to significance calculations using the MI- 
FAR statistic (see Appendix D). 


Appendix B: Event- by-event vetoes 

Often, GW candidate events identified in the on-source 
time can be linked to disturbances propagating through 
the detector from the environment or within the detector 
itself. Our procedure for identifying such event-by- event 
vetoes in S6 and VSR2/3 follows that used in S5 and 
VSR1 [16, 18]. The GW channel and a large number of 
auxiliary channels are processed with the Kleine- Welle 
[43] algorithm, which looks for excess power transients. 
A hierarchical method [44] is used to rank the statistical 
relationship between the transients found in the auxiliary 
channels and those found in the detector output. Based 
on these rankings, vetoes are defined for suspected noise 
events. Another veto used was based on significant sta- 
tistical association of events observed in the GW channel 
and the auxiliary channels [45] . 

An additional set of Category 3 vetoes [46] are applied 
to events emerging from networks including Virgo; vetoes 
from this set are based on detector read-out channels 
which are known to be insensitive to gravitational waves. 

Procedurally, the event-by-event vetoes are applied 
with the same conditions as their corresponding Cate- 
gory of data quality flags described in Appendix A. 


Appendix C: Detector networks and live times 

The total observation time for the analysis has been 
divided into four sub-periods (labeled A, B, C and D), 
separated by planned commissioning and upgrade breaks 
which changed the performance of the detectors. Ta- 
ble VI shows the observation time of each network config- 
uration after the application of Category 1 and 2 DQFs. 
These times are not overlapping. During the period from 
January to June 2010 (sub-period C), Virgo did not par- 
ticipate in the run because of hardware upgrades. 


detectors 

A [days] 

B [days] 

C [days] 

D [days] 

TOT [days] 

H1L1V1 

10.6 

16.7 

- 

24.9 

52.2 

H1L1 

- 

6.2 

51.4 

26.8 

84.5 

L1V1 

10.2 

10.7 

- 

8.1 

28.9 

H1V1 

12.6 

21.3 

- 

7.1 

41.0 

TOT 

33.4 

54.8 

51.4 

66.9 

206.6 


TABLE VI: Observation time for each detector configuration 
after application of Category 1 and 2 DQFs for the four sub- 
periods A, B, C, and D. For period A, the observation time of 
the H1L1 network after subtracting the H1L1V1 observation 
time is negligible (~ 1 day). During period C, Virgo did not 
participate in the run. 
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Appendix D: Modified Inverse False Alarm Rate 
(MIFAR) 

We introduce the Modified Inverse False Alarm Rate 
(MIFAR) to account for the effect of Category 3 DQFs 
on the background. 

Category 3 DQFs indicate a weak statistical correla- 
tion of the GW data with environmental and instrumen- 
tal noise sources, and thus were used only as a caution- 
ary tag when examining an event in candidate follow-ups. 
Moreover, the effectiveness of these flags is not constant 
between different sub-periods, network configurations or 
frequency bands. The use of Category 3 data quality as 
a tag allows us to produce two sets of events: the “raw” 
set (polluted to some extent by noise glitches) and the 
subse: of those events that are not tagged, the “clean” 
set (with reduced observation time). 

In order to account for the difference in background 
distributions when assessing the significance of candidate 
GW events from the raw and clean sets, we use the fol- 
lowing procedure: 

1 . Within each homogeneous analysis (same detec- 
tor’s configuration, same tuning of analysis, same 
frequency band), we rank events from the two sets 
separately by their coherent network amplitude rj ; 
i.e if the event candidate is flagged by Category 
3 data quality, it is ranked against the raw set of 
events, otherwise it is ranked against the clean set. 

2. Each event is then assigned a MIFAR as the inverse 


of the rate of higher-ranked background events in 
that set, i.e. the MIFAR is the IFAR of the event 
considering only that set. 

3. We merge the events from the raw and clean sets 
into a single list, sorted by the MIFAR. For events 
with equal MIFAR the one with larger 17 is ranked 
higher. This ranking is performed separately for 
each homogeneous analysis. 

4. According to this merged ranking, we measure the 
IFAR of the events as the rate of the corresponding 
background event with equal MIFAR. This mea- 
sured IFAR is used as our “universal ’ ranking for 
all events in all analvses. 

5. The final IFAR of any event over the entire search 
is just 1/3 of the value estimated in the previous 
step because of the trials factor: three independent 
analyses have been performed for the three disjoint 
frequency bands. No contribution to the trials fac- 
tor comes from the analyses of different detectors 1 
configuration since these were performed on non- 
overlapped observation times. 

In each homogeneous analysis, setting a threshold on 
IFAR corresponds to two thresholds on r) : one for the raw 
set and one for the clean data set. Table VII reports the 
selected thresholds. These thresholds were tuned using 
the background and injection events, without considering 
the on-source events to avoid bias in candidate selection. 
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network 

Free 

64-200 
cc 771 772 

uency Band [H 
200-2000 
CC 771 7)2 

z] 

2000-5000 
cc 771 772 

A 

H1L1V1 B 
D 

0.70 5.9 - 

0.60 6.5 6.5 
0.60 6.7 4.7 

0.70 6.0 - 

0.60 4.6 3.9 
0.60 4.5 4.3 

0.65 4.6 - 
0.60 3.9 - 
0.60 4.7 4.4 

A 

B 

H1L1 c 
D 

0.65 32.0 - 
0.60 8.5 8.1 
0.60 8.8 6.4 

0.60 8.9 8.9 

0.65 7.4 - 
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0.65 5.8 - 
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A 

H1V1 B 
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0.60 5.2 5.0 
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0.65 5.4 - 
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TABLE VII: Thresholds per network, sub-period and frequency band for the homogeneous analyses performed. The threshold 
rji is applied to the events in the raw set (those in coincidence with a Category 3 DQF) and 772 is applied to events in the clean 
set (net in coincidence with a DQF). These thresholds have been selected in order to ensure a IFAR > 8 yr. 
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